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c§ . ABSTRACT 



Recently the first detections of highly ionised gas associated with two Warm- Hot Intergalactic 
Medium (WHIM) filaments have been reported. The evidence is based on X-ray absorption lines 
«2j ■ due to O VII and other ions observed by Chandra towards the bright blazar Mrk 421. We 

d I investigate the robustness of this detection by a re-analysis of the original Chandra LETGS 

spectra, the analysis of a large set of XMM-Newton RGS spectra of Mrk 421, and additional 
Chandra observations. We address the reliability of individual spectral features belonging to the 
absorption components, and assess the significance of the detection of these components. We also 
use Monte Carlo simulations of spectra. We confirm the apparent strength of several features in 
the Chandra spectra, but demonstrate that they are statistically not significant. This decreased 
significance is due to the number of redshift trials that are made and that are not taken into 
account in the original discovery paper. Therefore these features must be attributed to statistical 
fluctuations. This is confirmed by the RGS spectra, which have a higher signal to noise ratio than 
the Chandra spectra, but do not show features at the same wavelengths. Finally, we show that 
the possible association with a Lya absorption system also lacks sufficient statistical evidence. 
We conclude that there is insufficient observational proof for the existence of the two proposed 
WHIM filaments towards Mrk 421, the brightest X-ray blazar on the sky. Therefore, the highly 
ionised component of the WHIM still remains to be discovered. 

Subject headings: BL Lacertae objects: individual: Mrk 421 - quasars: absorption lines - large-scale 
structure of Universe - X-rays: ISM 



1. Introduction 

Cosmological simulations show that the matter 
in the Universe is not uniformly distributed but 
forms a web-like structure. Clusters of galaxies 
and superclusters are found at the knots of this 
cosmic web. They are the location of the high- 
est mass concentration. The knots are connected 
through diffuse filaments, containing a mixture of 
dark matter, galaxies and gas. According to these 
simulations, the major part of this gas should be 
in a low density, intermediate temperature (10^ 
- 10^ K) phase, the so-called Warm-Hot Inter- 
galactic Medium (WHIM, Cen & Ostriker (1999)). 
About half of all baryons in the present day Uni- 
verse should reside in this WHIM according to the 
theoretical predictions, yet little of it has been seen 
so far. As the physics of the WHIM is complicated 
and many processes are relevant to its properties, 
observations of the WHIM are urgently needed. 

The density of the WHIM is low (typically 1 - 
1000 times the average baryon density of the Uni- 
verse (for example Dave et al. (2001)) or 0.3 - 
300 m^'^) and the temperature is such that most 
emission will occur in the EUV band. As the EUV 
band is strongly absorbed by the neutral hydro- 
gen of our Galaxy and its thermal emission is pro- 
portional to density squared, only the hottest and 
most dense part of the WHIM can be observed in 
emission with X-ray observatories such as XMM- 
Newton. Indications of O VII line emission from 
the WHIM near clusters have been reported now 
(Kaastra et al. (2003), Finogucnov ct al. (2003)). 

It is also possible to observe the WHIM in ab- 
sorption, provided that a strong background con- 
tinuum source is present and provided a high- 
resolution spectrograph is used. The first unam- 
biguous detections of the WHIM in absorption 
have been made in the UV band using bright 
quasar absorption lines from O VI as observed 
with the Far Ultraviolet Spectroscopic Explorer 
(FUSE) satelhte (Tripp et al. (2000); Oegerle et 
al. (2000); Savage et al. (2002); Jenkins ct al. 
(2003)). 

Observation in the UV band is relatively easy 
because of the high spectral resolution of FUSE 
and the relatively long wavelength A of the UV 
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lines. Since the optical depth at line centre scales 
proportional to the wavelength and the best spec- 
tral resolution of the current X-ray observatories 
is an order of magnitude poorer than FUSE, X-ray 
absorption lines are much harder to detect. 

The first solid detection of X-ray absorption 
lines at z = towards a quasar has been made by 
Nicastro et al. (2002). A spectrum of the bright 
quasar PKS 2155—304 taken with the Low Energy 
Transmission Grating Spectrometer (LETGS) of 
Ghandra showed zero redshift absorption lines due 
to O VII, O VIII and Ne IX. These lines have 
been confirmed by XMM-Newton Reflection Grat- 
ing Spectrometer (RGS) observations (Paerels et 
al. (2003); Cagnoni et al. (2003)), and similar X- 
ray absorption lines have been seen towards sev- 
eral other sources, for example 3C 273, Fang et 
al. (2003); NGG 5548, Stcenbrugge et al. (2003); 
NGC 4593, McKernan et al. (2003); and Mrk 279, 
Kaastra et al. (2004). However, whether these 
lines are due to the local WHIM (for instance 
around our Local Group) or that these lines have 
a Galactic origin is still debated. 

Therefore, the first unambiguous proof of X- 
ray absorption lines from the WHIM has to come 
from z > absorption lines. A first report of 
a z = 0.055 X-ray absorption system towards 
PKS 2155—304 has been given by Fang et al. 

(2002) using a Ghandra LETG/AGIS observation. 
However, this feature must be either a transient 
feature and therefore not associated to the WHIM, 
or an instrumental artifact. This is because these 
results have not been confirmed in both the more 
sensitive RGS observations towards this source 
(Paerels et al. (2003); Gagnoni et al. (2003)), 
nor in a long exposure using the LETG/HRG-S 
configuration. The X-ray detection of 6 interven- 
ing UV absorption systems towards H 1821+643 
in a 500 ks Ghandra observation by Mathur et al. 

(2003) should best be regarded as an upper limit 
due to the rather low significance of these detec- 
tions. 

An apparently more robust detection of z > 
X-ray absorption due to the WHIM has been pre- 
sented by Nicastro ct al. ((2005a), (2005b), NME 
herafter). These authors observed the brightest 
blazar on the sky, Mrk 421, when the source was in 
outburst, typically an order of magnitude brighter 
than the average state of this source. They ob- 
tained two high-quality Ghandra LETGS spectra 



of Mrk 421. They report the detection of two in- 
tervening X-ray absorption systems towards this 
source in the combined spectra. The first one at 
z = 0.011 with a significance of 5.8cr was associ- 
ated with a known H I Lya system, and the second 
one at z — 0.027 with a significance of 8.9(7 was 
associated with an intervening filament '^ 13 Mpc 
from the blazar. 

The cosmological implications of the detection 
of the first X-ray forest are perhaps the most im- 
portant clue to resolving the problem of the " miss- 
ing baryons" (see Nicastro et al. (2005a)), and 
therefore the robustness of this detection is ex- 
tremely important. Mrk 421 has been observed 
many times by XMM-Newton as part of its rou- 
tine calibration plan. The total integration time 
of 950 ks accumulated over all individual obser- 
vations with a broad range of flux levels yields a 
spectrum with a quality superior to the Chandra 
spectra taken during outburst. With such a high 
signal to noise ratio of the spectra all kinds of sys- 
tematic effects become important and we have de- 
veloped a novel way to analyse the RGS spectra 
taking account of these effects. This new analy- 
sis is presented in the accompanying paper (Ras- 
mussen et al. (2006)). Here we present a carcfuU 
re-analysis of all archival Chandra spectra of this 
source and compare the results from both instru- 
ments. 

2. Observations and data analysis 
2.1. Analysis of the LETGS spectra 

We processed all 7 Chandra LETGS observa- 
tions of Mrk 421 with a total exposure time of 
450 ks (see Table 1). Two deep (~ 100 ks) Chan- 
dra LETG/ACIS and LETG/HRC-S observations 
of Mrk 421 were triggered after outbursts catch- 
ing the source at its historical maximum. The 
^ 20 ks LETG/HRC-S observation was obtained 
in Chandra Director Discretionary Time after 2 
weeks of intense activity. The remaining four 
LETG/ACIS observations were performed for cal- 
ibration purposes. Our processed spectra contain 
in total ~ 10^ counts. There are also two short 
observations with the High Energy Transmission 
Grating Spectrometer (HETGS) of Chandra, but 
at the relevant wavelengths (A > 20 A) the signal 
to noise ratio of these spectra is too low for our 
present investigation. 



The LETG/ACIS spectra were processed with 
CIAO 3.2.2 (CALDB 2.27) using the standard 
Chandra X-ray Center pipeline. The spectra 
and ARFs were combined using the CIAO tool 
add-grating spectra. The LETG/HRC-S spectra 
were processed as described in detail in Kaastra 
et al.(2002). 

We combined these seven spectra into three 
datasets as indicated in Table 1: 

1. the ~ 100 ks LETG/ACIS observation per- 
formed in 2000 when the source had an ex- 
ceptionally high luminosity (the spectrum 
contains '^ 4.3 million photons) 

2. the combined two datasets obtained with 
LETG/HRC-S in 2000 and 2003 (the spec- 
trum contains ^ 2.7 million photons) 

3. the combined dataset of 4 calibration ob- 
servations with LETG/ACIS taken in 2004, 
when the source was on average at a lower 
luminosity (the spectrum contains ^3.2 mil- 
lion counts). 

We first analysed spectra 1 & 2, which are the 
data sets analysed by NME. However, contrary to 
NME we did not combine these data obtained in 
different instrument configurations (LETG/ACIS 
and LETG/HRC-S) and rather decided to fit the 
spectra simultaneously. For the spectral fitting we 
used the SPEX package (Kaastra et al. (1996)). 
The two spectra obtained with different instru- 
ment configurations were fitted with different con- 
tinuum models but with a common rcdshift and 
common absorption components. Since we are in- 
terested in weak absorption features and not in 
the continuum emission of the source, we fitted 
the underlying continuum with a spline model, us- 
ing 12 grid points between 9.5-35 A, with which 
we also removed any remaining broad instrumen- 
tal features (for more details on the spline model 
see the SPEX manual^). This spectral range for 
fitting was chosen because it includes all possible 
relevant spectral line features from neon, oxygen, 
nitrogen and carbon. 

We modelled the Galactic absorption using the 
hot model of the SPEX package, which calculates 
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Table 1: Chandra LETGS observations of Mrk 421. The cohinin "Set" gives the number of the spectral data 
set into which each observation was combined. The column "Flux" gives the flux at 20 A in phot m^^ s^^ A^^. 



Set 


Date 


Obs. 
ID 


Detector 


Exposure 

(ks) 


Flux 


1 


2000 May 29 


1715 


HRC-S 


19.84 


93 


2 


2002 Oct 26 


4148 


ACIS-S 


96.84 


333 


1 


2003 Jul 01 


4149 


HRC-S 


99.98 


242 


3 


2004 May 06 


5318 


ACIS-S 


30.16 


256 


3 


2004 Jul 12 


5331 


ACIS-S 


69.50 


39 


3 


2004 Jul 13 


5171 


ACIS-S 


67.15 


163 


3 


2004 Jul 14 


5332 


ACIS-S 


67.06 


140 



the transmission of a plasma in coUisional ioni- 
sation equilibrium with cosmic abundances. The 
transmission of a neutral plasma was mimicked 
by putting its temperature to 0.5 eV. We fix the 
Galactic absorption to 1.61 x 10^'* m~^ (Lockman 
& Savage (1995)). We fit the absorption lines of 
the z = local warm absorber, which is associated 
with the ISM of our Galaxy or the Local Group 
WHIM filament, also with the hot model of the 
SPEX package, but now with the temperature as 
a free parameter. Our best fit temperature of the 
local warm absorber is 0.07 keV with a hydrogen 
column density of 7.4 x 10^^ m~^. 
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Fig. 1. — Improvement Ax^ of the fit with weak 
absorption lines in the model with respect to a 
fit without these absorption lines, as a function 
of redshift. Ions included are C VI, N VI, N VII, 
O VI, O VII, Ne IX, and Ne X. The solid and 
dashed lines represent the Ax^ distribution as a 
function of redshift when fitting the observed spec- 
trum and the simulated spectrum respectively. 



After obtaining a good fit of the continuum 
and the local warm absorber we searched for the 
presence of weak absorption lines in the spectrum. 
We searched for redshifted absorption lines from 
C VI, N VI, N VII, O VI, O VII, O VIII, Ne IX, 
Ne X in all possible WHIM filaments between our 
Galaxy and Mrk 421. We performed a systematic 
search for these absorption features. Absorption 
lines of the above mentioned ions were placed into 
one of 60 bins corresponding to different redshifts 
between z = and z = 0.03. In each bin the in- 
dividual column densities of the ions were fitted 
and the Ax^ of the fit relative to the fit without 
the weak absorption lines in the model was eval- 
uated. To fit the column densities of the ions we 
used the slab model of the SPEX package. That 
model calculates the transmission of a thin slab 
of matter with arbitrary ionic composition. Free 
parameters are the velocity broadening av which 
we kept fixed to 100 kms~^, and the ionic column 
densities. Note that in our case the column densi- 
ties are low so we expect to be in the linear part of 
the curve of growth, such that CTv is irrelevant for 
the determination of the equivalent width of the 
absorption lines. In order to really detect weak 
absorption lines we performed this search with a 
spacing of 10 mA, one fifth of the resolution of 
the LETGS (50 mA). Afterward we simulated a 
Chandra spectrum with the same exposure time, 
same continuum model but without the weak ab- 
sorption features and again performed the same 
search on the simulated spectrum. The results of 
the search on the real and on the simulated spec- 
trum are shown in Fig. 1. 

We did not include redshifted O VIII fines, 
since a part of the spectrum where they may be 
present is influenced by small instrumental uncer- 



tainties. In particular a nearby node boundary in 
LETG/ACIS with small-scale effective area uncer- 
tainty would otherwise systematically bias the fit 
(see Sect. 3.1.1 for more details). 




o ^ 

o 



+ ++ 



^^■^vV^i 



i V^t^ l'l ^ l ^ l /tt/H l tV At^^^^ ^ l >t l W l t\ l tt^ l tVt l| t l t,h ^ 




I I H. I 4+ I +-'-* ! I l ^fi i 



21.5 22 

Wavelength (A) 



22.5 



Fig. 2.— The 21.0 - 22.6 A part of the spec- 
trum of Mrk 421 with its best-fitting continuum 
model plus an absorption model for the z = ab- 
sorber (solid line). The spectrum at the top is 
the LETG/ACIS observation obtained when the 
source was in an exceptionally bright state; the 
spectrum in the middle is the LETG/HRC-S ob- 
servation obtained in a bright state combined with 
a short LETG/HRC-S observation obtained at a 
lower luminosity state; these observations were 
combined and analysed by NME. The spectrum 
at the bottom was obtained during 4 calibration 
observations. The two LETG/ACIS spectra were 
multiplied by a factor of 1.5 for clarity of dis- 
play. The dotted line shows the position of the 
redshifted O VII line observed by NME. 

The most significant individual spectral feature 
which NME interpret as an absorption line from a 
WHIM filament is the feature associated with an 
O VII absorption line at z=0.011. To investigate 
that line further, wc considered all three datasets 
1-3. We show the 21.0 - 22.6 A spectral interval 
extracted from these datasets in Fig. 2. The best- 
fitting continuum model plus an absorption model 
for the local warm absorber is shown with a solid 
line and the position of the redshifted O VII line 
observed by NME is indicated by the vertical dot- 
ted line. The spectral feature at 21.85 A is seen 
in the spectrum obtained by LETG/HRC-S and 
a somewhat weaker and shifted feature is seen in 



the LETG/ACIS spectrum obtained at the high 
source luminosity. We note that NME combined 
these two datasets. In the third dataset we do 
not see any feature at the indicated wavelength, 
despite the fact that the spectrum contains more 
photons than the LETG/HRC-S observation. 

Finally, we determined the nominal rcdshift of 
the z = 0.011 component from the data sets 1 and 
2, using only the strongest absorption lines (O VII 
ls-2p). For the z = O VII line, we measure a line 
centroid of 21.605 ± 0.006 and 21.604 ± 0.005 mA 
for spectrum 1 and 2, respectively; for the z — 
0.011 components these values are 21.864 ± 0.014 
and 21.854 ± 0.011 mA. Therefore, the redshifts 
of the longer wavelength line with respect to the 
z = line is 0.0121 ± 0.0007 and 0.0117 ± 0.0005, 
respectively. The weighted average is z = 0.0118± 
0.0004, somewhat higher than the value obtained 
by NME. 

2.2. Analysis of the RGS spectra 

The analysis of the RGS spectra is described 
in detail by Rasmussen et al. (2006). We have 
taken their combined fluxed spectra and corrected 
these for Galactic absorption using the same hot 
model of SPEX described above that was used for 
the LETGS spectra. We adjusted the oxygen col- 
umn density to 8.5 x lO"** of the hydrogen column 
density in order to get the best match around the 
oxygen edge. The absorption corrected spectra 
then were fitted using a spline model with nodes 
separated by 0.5 A. Obvious absorption lines were 
ignored in this continuum fit. This procedure is 
needed because the statistical quality of our com- 
bined data is so good: typically, the flux in the 
oxygen region is determined with a statistical er- 
ror of only 0.4 % per 0.5 A bin. There are re- 
maining, large-scale uncertainties in the RGS ef- 
fective area of the order of a percent on A scales 
that would otherwise bias the spectrum, and more- 
over to prove that the underlying continuum of all 
these combined Mrk 421 observations would be a 
straight power law would be a major challenge. 
Fig. 3 shows the residuals of this fit in the same 6 
wavelength intervals used by NME. 

In addition, we show in Fig. 3 the fit residuals 
of the LETGS spectra 1 and 2, after fitting with 
the same spline model with with nodes separated 
by 0.5 A as was used for the RGS data. 
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Fig. 3.— Comparison of the LETGS and RGS spectra of Mrk 421 in the Ne X/Ne IX, O VIII, 
O VII/0 VI, N VII, N VI and C VI bands. Each panel shows from top to bottom the LETG/ACIS 
spectrum (1), the LETG/HRC-S spectrum, RGSl and RGS2. The quantities plotted are the fit residuals 
(observed— model) /model with respect to a local continuum model. One tickmark on the y-axis corresponds 
to 10 %, and the data sets have been shifted arbitrarily along the y-axis in order to avoid overlapping data 
points. The solid vertical lines correspond to expected features at 2; = 0, dashed lines to z = 0.011, and 
dash-dotted lines to 2; = 0.027. CCD gaps for RGS are shown with dotted lines in the appropriate spectrum. 
As explained by Rasmussen et al. (2006), after careful screening occasionally an isolated instrumental dip 
may be present in a single RGS, such as at 25.8 A in RGS2. 



2.3. Determing equivalent widths 

Equivalent widths of selected line features in 
the 21.0-22.5 A band were determined from the fit 



residuals shown in Fig. 3. We consider the z = 
line of O VII, as well as the z — and z = 0.027 
O VII lines and a 22.02 A feature discussed later. 



Table 2: The formal best-fit equivalent widths (W\) in niA for absorption lines fixed at wavelengths of 21.60, 
21.85, 22.20 and 22.02 A, corresponding to the expected position of the ls-2p lines of O VII at z = 0, 
z = 0.011, z = 0.027, and the ls-2p doublet of O VI at z = 0. Negative values indicate emission lines. See 
the text for the details about the LETGS datasets. We further list W\ as derived from the RGS spectra 
(Sect. 2.2), and the weighted average for all three LETGS data sets and the RGS data (row labelled w.a.). 
For comparison, we also list the equivalent widths as determined by Williams et al. (2005) and NME from 
the combined LETGS sets 1 and 2 (labelled NME/W). All errors are r.m.s. errors (Ax^ = 2). The last 
column gives an estimate of the systematic uncertainty (also in niA) on the equivalent widths (not included 
in the statistical errors). 



Datasct 


21.60 


21.85 


22.20 


22.02 


syst 


LETGS 1 
LETGS 2 
LETGS 3 
RGS 


11.6 ±1.3 
11.7±1.3 
10.3 ±1.6 
14.8 ±0.7 


1.5 ±1.5 

3.5 ±1.4 

-0.6 ±1.8 

0.4 ±0.8 


3.1 ±1.5 
-1.5±1.5 
-0.4 ±1.9 
-0.6 ±0.7 


1.9 ±1.5 
1.1 ±1.5 
0.1 ±1.8 
2.8 ±0.7 


0.6 
0.9 
0.6 

0.8 


w.a. 


13.4 ±0.5 


1.0 ±0.6 


-0.2 ±0.6 


2.1 ±0.6 




NME/W 


9.4 ±1.1 


3.0 ±0.9 


2.2 ±0.8 


2.4 ±0.9 





The equivalent widths {W\) are given in Table 2. 

For RGSl we used the exact line spread func- 
tion. Following our earlier work (Kaastra et al. 
(2002)), we approximate the Isf of the LETGS by 
the following analytical formula: 



(f>^l/[l + (AX/af 



(1) 



where a = 2 and AA is the offset from the line 
centroid. The scale parameter a (expressed here 
in A) is weakly wavelength dependent as: 



a= \/ (0.0313)2 + (0.000306A)2 



(2) 



The errors on the equivalent widths that we list 
are simply the statistical errors on the best-fit line 
normalisation. Table 2 also shows the systematic 
uncertainties expected for each line due to the un- 
certainty in the continuum level; these systematic 
uncertainties were determined by shifting the con- 
tinuum up and down by the la statistical uncer- 
tainty in the average continuum flux for each local 
0.5 A bin and re-determing the equivalent width. 

Uncertainties in the shape of the Isf also may 
affect the derived values of W\. We estimated 
this by modifying the Isf and determining how the 
equivalent width of the 21.60 A line changes. For 
the LETGS, a change of a from 2 (our preferred 
value) to 2.5 (the value used by the official CIAO 
software) but keeping the same FWHM by adjust- 
ing a leads to a decrease of 0.4 mA (3.5 %)in Wy, a 
decrease of a by 10 % leads to a decrease of 0.5 mA 



(4.5 %)in Wx. If for the RGS instead of the true 
Isf a Gaussian approximation would be used, the 
equivalent widths would be 20-30 % too small. As 
shown by Rasmussen et al. (2006) the uncertainty 
in the Isf gives slightly different results (^ 10 %) 
as two different parameterisations of the Isf are 
used. For example, decreasing the FWHM of the 
Isf of the RGS by 10 % leads to a W\ decrease of 
2.3 mA (15 %). Although the precise systematic 
uncertainty on W\ due to the Isf shape is hard to 
establish accurately, it is most likely in the 5-15 % 
range for both LETGS and RGS. Taking this into 
account in addition to the statistical uncertainty 
on Wx and the systematic uncertainty related to 
the continuum level (Table 2), we conclude that 
there is not a large discrepancy between the equiv- 
alent widths of the z = O VII line measured by 
both instruments. 

3. Discussion 

3.1. Significance of the absorption compo- 
nents detected by Chandra 

The evidence for the two absorption compo- 
nents at z = 0.011 (5.8cr) and z = 0.027 (8.9ct) 
as presented by NME looks at a first glance con- 
vincing. However, there are good reasons to put 
question marks to these significances. We follow 
here two routes: first we discuss the significance 
by looking into detail to the individual features 
identified by NME, and next we analyse the signif- 



icance of the absorption components in a coherent 
way. 

3.1.1. Individual line features 

NME found 28 or 24 absorption features with 
a significance of more than 3cr in their analysis 
of the combined LETG/HRC-S and LETG/ACIS 
spectrum. The precise number is somewhat un- 
clear from their paper but is not important here. 
Of these 24 lines, 14 have an origin at z == and 
there is no doubt about the significance of that 
component, so we focus on the 10 features with 
a tentative identification with z > absorption 
lines. In total 3 of these belong to the z — 0.011 
system, and 6 to the z — 0.027 system. One fea- 
ture at 24.97 A was unidentified by NME. We will 
discuss these features here individually as we point 
out that not all of these detections can be regarded 
as unambiguous. 

The absorption component at z = 0.011 had 
three features: 

1. O VIII ls-2p at 19.18 A: Fig. 3 shows that 
there is only evidence for an absorption 
feature from the LETG/ACIS spectrum, 
not from the LETG/HRC-S or RGS spec- 
tra. The feature coincides exactly with a ~ 
0.05 A wide dip of 3 % depth in the effective 
area of the LETG/ACIS, associated with an 
ACIS node boundary (as was also pointed 
out by NME); therefore this feature is most 
likely of instrumental origin. 

2. O VII ls-2p at 21.85 A: there is no clear fea- 
ture in the RGS spectrum, but there is a fea- 
ture in both LETGS spectra. However, the 
centroid in the LETG/ACIS spectrum seems 
to occur at slightly longer wavelength, by 
~ 0.01 A. There are no known instrumental 
features in the LETGS at this wavelength, 
so this may be real provided it has sufficient 
significance. But see Sect. 3.1.2. 

3. N VII ls-2p at 25.04 A: no feature present in 
the RGS deeper than 1 %; the LETG/HRC- 
S spectrum (Fig. 3) shows a '--^ 5 % negative 
deviation at a slightly smaller wavelength, 
and the (here noisy) LETG/ACIS spectrum 
may be consistent with this. The smaller 
wavelength is consistent with the findings of 



NME, who give here a redshift of 0.010 in- 
stead of 0.011. 

In summary, the X-ray evidence for the 0.011 
component rests solely upon the detection of two 
lines with different redshift and only seen in the 
LETGS. As we will argue later, even the statistical 
significance of the combined features is insufficient 
to designate it as a robust detection. 

The component at z == 0.027 had six features: 

1. Ne IX ls-2p at 13.80 A: this feature coin- 
cides with a narrow dip in the effective area 
of the LETG/ACIS associated with a node 
boundary, similar to the z ~ 0.011 O VIII 
ls-2p line. It is therefore no surprise that it 
is only seen in the LETG/ACIS spectrum, 
with perhaps a little but not very significant 
strengthening from the LETG/HRC-S spec- 
trum. It is not visible in the RGS spectrum. 
More importantly, its best-fit redshift as re- 
ported by NME of 0.026 is off from the red- 
shift of this component at z = 0.027. We 
conclude that this feature has an instrumen- 
tal origin. 

2. O VII ls-3p at 19.11 A: there is only a weak 
signal visible in the LETG/HRC-S spectrum 
(Fig. 3); this looks a little less pronounced 
than in the plot given by NME but this may 
be due to different binning; more impor- 
tantly, given the weakness of the ls-2p tran- 
sition in the same ion (see below) and the 
expected 5 times smaller equivalent width of 
the ls-3p line as compared to the ls-2p line 
makes any detection of the ls-3p line impos- 
sible. 

3. O VII ls-2p at 22.20 A: only clearly seen 
in the LETG/ACIS spectrum. No evident 
instrumental features are present. However, 
the redshift of 0.028 is slightly off from the 
redshift of this component at z == 0.027. 

4. N VII ls-2p at 25.44 A: a weak feature is 
present mainly in the LETG/HRC-S spec- 
trum; RGSl may be consistent with this. 
Given the weakness of the oxygen lines, a 
clear detection of N VII would be surprising 
unless the absorber has an anomalously high 
nitrogen abundance. 



5. N VI ls-2p at 29.54 A: here the LETG/HRC- 
S spectrum has the best statistics and it 
shows indeed a shaUow dip with four adja- 
cent bins approximately la befow the ref- 
erence level. The same remark as above 
about the nitrogen abundance can be made. 
Moreover, there is nothing there in the RGS 
spectra. 

6. C VI ls-2p at 34.69 A: as above the feature is 
mainly evident from the LETG/HRC-S line. 
As the cosmic carbon abundance is signifi- 
cantly higher than the nitrogen abundance, 
this line is more likely to be detectable. 
Again, the redshift of 0.028 is on the high 
side. 

In summary, the detection of the z ~ 0.027 
component rests upon O VII, seen only in the 
LETG/ACIS spectrum, and the ls-2p lines of 
N VII, N VI and C VI as well as the ls-3p lines of 
O VII, seen only in the LETG/HRC-S spectrum. 
None of these features is visible in both LETGS 
configurations, not all of these lines have exactly 
the same redshift and some of the lines (O VII Is- 
3p and the nitrogen lines) are problematic from a 
physical point of view. 

3.1.2. Self- consistent assessment of the signifi- 
cance of the line features 

In the previous section we have put several 
question marks to the significance and reliability 
of the nine z = 0.011 and z = 0.027 line identi- 
fications as derived from the LETGS spectra by 
NME. Nevertheless, some of the features identi- 
fied by these authors indeed show negative resid- 
uals at the predicted wavelengths for these com- 
ponents, with a nominal statistical significance as 
determined correctly by NME. 

However, the actual significance of these detec- 
tions is much smaller, and actually as we show 
below the redshift components are not significant 
at all. This is due to two reasons. 

First, NME fitted line centroids for each com- 
ponent separately, and then determined the sig- 
nificance for each line individually. However, for 
all lines belonging to the same redshift compo- 
nent, the wavelengths are not free parameters but 
are linked through the value of the redshift. Al- 
though NME later assess the significance of the 



components by fixing the wavelengths of the rel- 
evant lines, they first have established the pres- 
ence of redshifted components and their redshift 
values by ignoring this coupling. Their argument 
for "noise" in the wavelength solution of in par- 
ticular the LETG/HRC-S detector is not relevant 
here; first, our improved treatment of the LETGS 
gives smaller noise, but more importantly, if there 
is noise, then for each potential line the amplitude 
and direction of the noise are unknown, and there- 
fore there is no justification for adjusting the line 
centroids by eye or by chi-square: the true wave- 
length alinearity might just have the opposite sign. 
The only way to do this unbiased is to couple the 
lines from different ions even before having found 
the redshifted components. 

The most important problem is however the 
number of trial redshifts that have been consid- 
ered. NME started by looking to > 3(7 excesses 
in the data in order to identify the two redshift 
systems. However, a > 3cr excess has only mean- 
ing for a line with a well-known wavelength based 
on a priori knowledge. For example, knowing 
that there is a highly ionised system with a given 
redshift in the line of sight makes it justified to 
look for O VIII or O VII absorption lines at pre- 
cisely that redshift. Here this a priori knowledge is 
lacking (only a posteriori circumstantial evidence 
based on galaxy association or association H I Lya 
absorption is presented) . In a long enough stretch 
of data from any featureless spectrum there will 
always be several statistical fiuctuations that will 
surpass the 3ct significance level. The significance 
of any feature with unknown wavelength is there- 
fore not determined by the cumulative distribution 
F{x) of the excess x (with F{x) usually a Gaus- 
sian) at a given wavelength. Instead it is given by 
the distribution G{x) of the maximum of F for N 
independent trials, which is G{x) = F'^{x). For 
a large number of trials, F has to be very high to 
get a significant value of G. 

The number of independent trials A^ is of or- 
der of magnitude W/AX, with W the length of 
the wavelength stretch that is being searched, and 
AA the instrumental Full- Width at Half Maximum 
(FWHM). However, this is only an order of mag- 
nitude estimate. The precise value depends on the 
shape of the instrumental line spread function and 
the bin size of the data. 

By using a Monte Carlo simulation we have as- 



sessed the significance of line detections for the 
present context. We start from a simple feature- 
less model photon spectrum. This spectrum is 
binned with a bin size i5 of 0.015 A. For each spec- 
tral bin i, the deviations di of the observed spec- 
trum from this model spectrum are normalised to 
their nominal standard deviations Ci such that 
Xi = di/ai has a standard normal distribution 
with mean and variance 1. Given the high count 
rate of the spectrum, the approximation of the 
Poissonian distribution for each bin with a Gaus- 
sian as we do here is justified. All Xi are statis- 
tically independent random variates. We consider 
7 spectral lines, centred on \j equal to 12.1, 13.4, 
21.6, 22.0, 24.8, 28.8 and 33.7 A, and generate 
the Xi for the wavelength range between \j and 
Aj(l + Zmax) with Zmax the rcdshift of Mrk 421 
which we approximate here by 0.03. Then we do a 
grid search over redshift (step size 0.0003 in z) and 
for each redshift we determine the best fit absorp- 
tion line fluxes. The fiuxes are simply determined 
from a least squares fit taking account of the exact 
line profile 4>{X) (Eqn. 1). 

As we are looking for absorption lines, we put 
the normalisation of the line to zero whenever it 
becomes positive. We then calculate the difference 
A in x^ between the model including any absorp- 
tion lines found, to the model without absorption 
lines. As the fit always improves by adding lines, 
A% is always negative. For each run, we deter- 
mine the redshift with the most negative value of 
Ax^, i.e. the most significant signal for any red- 
shift. We repeat this process 10® times and de- 
termine the statistical distribution of Ax^ by this 
way (see Fig. 4). 

Our results are almost independent of the cho- 
sen bin size 5 for the spectrum (we verified that 
by numerical experiments). In the context of 
Mrk 421, the median and peak of the distribution 
are close to 10. This corresponds to a "5cr" de- 
tection if no account is taken of the number of 
trials. Components with apparent significances 
stronger than 5.8 or 8.9cr as obtained by NME 
for the z — 0.011 and z — 0.027 systems, respec- 
tively, have probabilities of 40 and 6 % to occur by 
chance. Thus, the z — 0.011 component is not sig- 
nificant at all. Furthermore, as we show below in 
the next section, the z = 0.027 is not as significant 
as reported by NME. 



3.1.3. The significance of the line spectra from 
spectral simulations 

The two peaks at z = 0.011 and z — 0.027 
are indeed reproduced by our new analysis of the 
LETGS spectra (Fig. 1; Sect. 2.1). While the peak 
aX, z = 0.011 has a significance similar to the value 
quoted by NME, the z — 0.027 component is less 
significant. This is mainly because we omitted the 
O VIII wavelength range because of the reasons 
mentioned earlier. Note that both troughs arc not 
very narrow; this corresponds to the slightly dif- 
ferent redshifts found for the lines of NME. 

We also took the same continuum spectrum 
and used it to make a simulation (Fig. 1). The 
simulated spectrum shows three peaks of similar 
strength, dX z = 0.001, z = 0.007, and z = 0.019. 
In both the fitted true spectrum of Mrk 421 and 
the simulated spectrum, the strongest peak in 
terms of A^^ is in excellent agreement with the 
predictions as shown in Fig. 4. 

We conclude that based on the Chandra data 
alone there is not sufficient evidence for a signif- 
icant detection of the WHIM towards Mrk 421. 
This does of course also imply that a true WHIM 
signal must be much stronger or more significant 
than the detections of NME before it can be re- 
garded as evidence for the existence of the WHIM. 

3.2. Association of the z = 0.011 compo- 
nent with a Lya absorber. 

A strong argument in favour of the existence 
of the z = 0.011 WHIM component is its associa- 
tion with a Lya absorption system. Penton et al. 
(2000) found an absorption line with an equivalent 
width of 86 ±15 mA at 1227.98 A, identified as H I 
Lya at a redshift of z = 0.01012 ±0.00002. Savage 
et al. (2005) have studied extensively UV spectra 
of Mrk 421 obtained by FUSE and HST but they 
could not find any other H I lines neither any as- 
sociated metal lines. However, Savage et al. argue 
that this line could be hardly anything else than 
Lya. At the given redshift, the system should be 
a H I cloud in a galactic void. 

From our analysis of the apparent z ^ 0.011 
system (Sect. 2.1), wc found as a best fit red- 
shift z = 0.0118 ± 0.0004. This is higher than the 
value quoted by NME because of our improved 
wavelength solution of the LETGS. The difference 
of 500 ± 120 kms~^ between this X-ray redshift 
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and the Lya redshift corresponds to 0.036 A at 
the wavelength of the O VII resonance hne. This 
is almost equal to the spectral resolution of the 
LETGS at that wavelength and definitely inconsis- 
tent with our current understanding of the LETGS 
wavelength scale. 

We conclude that there is no evidence for a 
direct physical connection between the putative 
z = 0.011 X-ray absorption component and the 
Lya absorber. A similar conclusion was also 
reached by Savage et al. (2005), based on the 
narrowness of the Lya line which implies a low 
temperature and the absence of any O VI absorp- 
tion in the FUSE spectra. Even the evidence for a 
more loose association of both systems, for exam- 
ple in different parts of a larger filamentary struc- 
ture, is not convincing: over the redshift range of 
< z < 0.03, the probability that the most sig- 
nificant statistical X-ray fluctuation would occur 
within \/S.z\ < 0.0017 from the only known inter- 
vening Lya system is 11.2 %. 

3.3. The final argument: no lines in the 
RGS spectrum 

We have shown above that based on the Chan- 
dra data alone there is insufficient proof for the 
existence of the z = 0.011 and z — 0.027 absorp- 
tion components reported by NME. The lack of 
evidence does not imply automatically that these 
components are not present, but just that their 
presence cannot be demonstrated using LETGS 
data. 

However, if these components really would ex- 
ist at the levels as reported by NME, then adding 
more data like our combined RGS data set should 
enhance the significance, but the opposite is true: 
the RGS spectra show no evidence at all for both 
absorption components. Taking the weighted av- 
erage of the equivalent widths of the O VII lines 
obtained by all instruments, we find for the z = 
0.011 and z — 0.027 components equivalent widths 
of 1.0 ± 0.6 and —0.2 ± 0.6 mA, respectively (see 
Table 2), to be compared to 3.0 ±0.9 and 2.2 ±0.8 
as reported by NME. Hence, the column densities 
of any component are at least two times smaller 
than reported by NME. 

The lack of evidence for both components in 
RGS spectra was also reported by Ravasio et al. 
(2005), based on a smaller subset of only 177 ks. 



Our analysis uses 4 times more exposure time. 

Finally, while we were writing our paper, a 
preprint by Williams et al. (2006) appeared. 
These authors use 437 ks RGS data, only half of 
our exposure time. They also do not find evidence 
for the two absorption systems, but attribute this 
to "narrow instrumental features, inferior spectral 
resolution and fixed pattern noise in the RGS" . As 
we have shown in the accompanying paper (Ras- 
mussen et al. (2006)) these statements lack justi- 
fication when a careful analysis is done. 

3.4. The 22.02 A feature 

Both the RGS and the LETGS spectra seem 
to indicate the presence of an absorption line 
at 22.02 A, with combined equivalent width of 
2.1 ± 0.6 mA (Table 2). This feature was identi- 
fied as z = absorption from O VI, based upon the 
almost exact coincidence with the predicted wave- 
length of that line (Wilhams et al. (2005)). How- 
ever, as Williams et al. point out, the equivalent 
width is about 3 times higher than the value pre- 
dicted from the 2s-2p doublet at 1032 and 1038 A. 
Williams et al. then offer three possible explana- 
tions for the discrepancy. 

The first explanation, wrong oscillator strengths 
of the K-shell transitions by a factor of 2-4, can 
be simply ruled out. These strong transitions can 
be calculated with much higher precision. 

The second explanation invokes a large fraction 
of O VI in an excited state. This would suppress 
the 2s-2p UV lines relative to the ls-2p X-ray 
lines. Apart from the fact that this situation is 
hard to achieve in a low density plasma (one needs 
almost LTE population ratio's), there is another 
important argument against this. The X-ray ab- 
sorption lines from the Is^ 2p configuration are 
shifted by 0.03-0.05 A towards longer wavelength 
as compared to lines from the ground state Is^ 2s 
(A.J.J. Raassen, private communication). Thus, 
in this scenario the X-ray line should shift by a 
measurable amount towards longer wavelengths, 
which is not observed. 

The third explanation offered by Williams et al. 
invokes intervening O VII a.t z — 0.0195. Given 
the lack of other absorption lines from the same 
system, Williams et al. argue that this is unlikely. 

Our strongest argument against a real line is 
however its significance. Using the known column 
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density derived from the FUSE spectra, (2.88 ± 
0.14) X 10^*^ ni-2, Savage et al. (2005)), we 
estimate that the corresponding z — ls-2p 
line of O VI should have an equivalent width 
of 0.64 mA. Subtracting this from the observed 
equivalent width, the "unexplained" part has an 
equivalent width of 1.5±0.6 mA, i.e. a 2.5ct result. 
However, similar to our analysis of the Chandra 
data, it is easy to argue that given the number of 
redshift trials this significance is insufficient to be 
evidence for redshifted O VII. 

4. Conclusions 

We have re-analysed carefully the Chandra 
LETGS spectra of Mrk 421, the first source for 
which the detection of the X-ray absorption for- 
est has been claimed. We have supplemented this 
with an analysis of additional LETGS observa- 
tions as well as with a large sample of XMM- 
Newton RGS observations, amounting to a total 
of 950 ks observation time. We find that there is 
no evidence for the presence of the z ~ 0.011 and 
z — 0.027 filaments reported by NME. Also the as- 
sociation with an intervening H I Lya absorption 
system is not sufficiently supported. 
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Fig. 4. — Probability density function (lower 
panel) and cumulative probability density func- 
tion (upper panel) for the maximum improvement 
— Ax^ in LETGS fits including absorption lines at 
12.1, 13.4, 21.6, 22.0, 24.8, 28.8 and 33.7 A, when 
a redshift range between z ~ and z — z^ax is 
searched. Here results for 3 values of z^ax are 
shown: dotted line: Zmax = 0.015; thick solid line: 
^max = 0.03 (appropriate for Mrk 421); dashed 



line: 



0.06. 
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